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We study translationally invariant rings of qubits with a finite number of sites N, and find the 
maximal nearest-neighbor entanglement for a fixed z component of the total spin. For small numbers 
of sites our results are analytical. The use of a linearized version of the concurrence allows us to 
relate the maximal concurrence to the ground state energy of an XXZ spin model, and to calculate 
it numerically for N < 24. We point out some interesting finite-size effects. Finally, we generalize 
our results beyond nearest neighbors. 
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INTRODUCTION 



Entanglement in many-body systems such as spin chains has recently attracted much attention. The reasons are at 
least twofold: on one side one is interested in understanding fundamental structures of entangled multipartite states 
and the natural occurrence of such states in systems with spin-spin interactions, with quantum information processing 
purposes in mind. On the other side insights from quantum information theory may lead to a deeper understanding 
of many-body phenomena in strongly interacting systems, such as quantum phase transitions and critical behavior. 
Much effort has been devoted to studying the natural appearance of bipartite entanglement in ground states of 1- 
| dimensional spin models with an Ising Hamiltonian or a Heisenberg-type Hamiltonian 0,07 with special emphasis 
00 ' on properties of entanglement near a quantum phase transition. It has been pointed out that some entanglement 
remains for finite temperatures (so-called thermal entanglement 5] ) , which for some transverse field can even increase 
with increasing temperature. Recently, it has been shown that maximally entangled states of atoms can be created 
by crossing a quantum phase transition 



A complementary approach to the phenomenon of entanglement in chains of qubits has been initiated by O'Connor 
and Wootters (see also @): they ask (without specifying a Hamiltonian) the fundamental question "what is 
the maximal entanglement between two neighboring sites of an entangled ring with translational invariance"? An 
"entangled ring" is defined as a chain of qubits with periodic boundary conditions. Due to the so-called "monogamy 
-y ' of entanglement" it is impossible for a site to be maximally entangled with both its neighbors: shared entanglement 
is always less than maximal 0, Il0| . In Ref. 7] the question of the upper limit for the nearest neighbor entanglement 
was simplified by introducing two additional restrictions on the allowed states: 



(i) the state is an eigenstate of the z component of the total spin, 

(ii) neighboring spins cannot both be "up" . 



Both restrictions are based on an educated guess for the optimal states for the general problem. O'Connor and 
Wootters solved the restricted optimization problem by relating it to an effective Hamiltonian for the one-dimensional 
ferromagnetic XY model, and found the maximal nearest-neighbor concurrence to be 

2sinf 1 ^- 

c max (A, P ) = — y_*<_ , (i) 

N sin I 



(iV-p) 



where N is the number of sites and p is the number of up-spins. Eq. provides a lower bound for the problem 
without restriction (ii). For given N and p, it may or may not happen that C can be increased by also allowing states 
where two neighboring spins are up. In the limit N — > oo, the optimal number of up-spins in Eq. Q approaches 
Popt ~ 0.301 N. This leads to an asymptotic value of C" nax w 0.434. Whether this number can be improved by 
removing restriction (ii) or restriction (i) is an open problem. We will illustrate, however, that for a fixed p restriction 
(ii) tends to play a decreasing role as N is increased. 
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In this article we mainly focus on the problem of optimizing the nearest-neighbor entanglement for given finite N 
and p, without imposing restriction (ii) . Using the methods described below we solve a number of cases analytically. 
We supplement these results by numerical calculations, and thus arrive at a rather complete picture of the solutions 
for N < 24. Some finite size effects that occur only for specific small numbers of sites will be pointed out. 

The paper is organized as follows: In Sec. |il|we introduce the problem and our notation, and study symmetry 
under reflections. In Sec. IHII we describe our analytical results for the maximal nearest-neighbor entanglement in the 
cases p = 2 and p = N/2. In Sec. IIVI we numerically investigate whether the solutions in are at least still local 
maxima when restriction (ii) is removed. In Sec. we use a method recently introduced by Wolf et al. mapping 
the original nonlinear problem onto a series of linear ones. We discuss some general properties of the mapping, and 
use it to numerically calculate the maximal nearest-neighbor concurrence for N < 24. Finally in Sec. IVUI we discuss 
our results and point out some interesting questions that are still open. 

II. MAXIMAL NEAREST NEIGHBOR ENTANGLEMENT 

We consider a ring of N qubits, out of which p are in the state | f) and N —p are in the state | J.). In the following 
we will use the notation | J.) = |0) and | f) = |1), and will refer to p as the occupation number. We will always 
consider translationally invariant states. 

A. Translationally invariant states and concurrence 

Our aim is to calculate the maximal nearest-neighbor entanglement of two qubits, using the concurrence as a 
measure. The concurrence is defined as C(p) = max{Ai — A 2 — A 3 — A4, 0}, where Ai > A2 > A3 > A 4 > are the 
square roots of the eigenvalues of pp, and p :— (a v <g> a v )p(a y (8 a v ) is the spin-flipped density matrix. 

The structure of any two-qubit reduced density matrix on the ring is, due to the properties of translational invariance 
and the occupation number p being fixed, given by 
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with v + 2w + y = 1. The basis is chosen as (|00), |01), |10), |11)). Here the zero entries are due to considering a fixed 
p. The diagonal entries in the density matrix are real, and the equality of two of the diagonal elements is caused by 
the translational invariance. Note that v — y + (N — 2p)/N, i.e. the number of down-down pairs of neighboring spins 
depends on the number of up-up pairs. We will see below that for states with maximal entanglement the off-diagonal 
element z can be taken to be real. 

The concurrence of the density matrix J2Jl is given by 

C = 2max{|z| - y/vy^} . (3) 

As one observes from Eq. the task of finding the state that maximizes the nearest- neighbor concurrence is 
substantially simplified by requiring y = 0, i.e. two neighboring spins cannot be both "up" . This is the assumption 
that was made in — we will relax this constraint throughout our work. 

For given N and p we apparently need to maximize Eq. J2J over all translationally invariant states, that is, over 
all density matrices fulfilling [g, T] = 0, where T is the operator that implements a translation by one site. As was 
shown in [jj, the problem splits into N maximization problems over pure states: since g and T commute, the spectral 
decomposition of g can be expressed in projectors onto eigenstates of T, namely 

JV-l 

Q=^Pk\il>k){^kl (4) 

fc=0 

where pk > and J^kPk — 1- Here the eigenstates of the translation operator are defined by T\ipk) = e l ^ k \ipk)- 
Since C is a convex function, we have 



C 5> fc |VkXVk| < 5>kCW*Wfc|) < max{C(|^)(^.|)}, 

V k / k 



(5) 
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and it is thus sufficient to first maximize the concurrence over pure states within each eigenspace of the translation 
operator, and then choose the maximal result over all eigenstates. In fact, the situation is even simpler, since, as 
shown in 0, the search can be restricted to the k = eigenstates, i.e. to eigenstates of T with eigenvalue 1. 

For a given number of sites JV and a given occupation number p, the most general eigenstate of T with eigenvalue 
1 is a superposition that looks as follows: 

\i>(N,P)) =E a » > 



with the normalizations 



and with 



(/i| M ) = l and ^2\a„\ 2 = l, (7) 



N-l 

W=^^T j W. (8) 

3=0 

The above notation is as follows: the index \i enumerates all distinct translationally invariant configurations of the 
ring where the p up-spins and (JV — p) down-spins have given relative positions. The state is one representative 
member of the /ith configuration. The operator inducing a translation by j sites is given by T- 7 , and Af^ denotes an 
appropriate normalization factor. Thus, is a normalized eigenstate of T corresponding to the configuration p, and 
the eigenvalue 1. Note that for a given /j, it may happen that T J |^>„) = already for j = JV < JV; in that case it 
is not possible to construct eigenstates corresponding to all eigenvalues of T from \4>n)- This is accounted for by an 

appropriate value of JV^, which in the general case (JV < JV) reads JV M = \/jV /JV. 

We use the following convention to denote /i explicitly: the composite index /j, consists of p entries, where each 
entry specifies a distance from one up-spin to the next up-spin in |</> M ). We choose the first entry to be the smallest. 
Note that the p entries add up to JV. An example will clarify this notation: let us consider the case of JV = 5 qubits 
and the occupation number p — 2. Then the index fi can only take two combinations of entries, namely 1,4 and 2,3. 
The most general translationally invariant state with translational eigenvalue 1 is written as 

|V(JV = 5,p = 2)) = oi,4|l,4) + o 2l 3|2,3) 

- E <^E ri i^> 

= ax 4^7= (|11000) + |01100) + |00110) + |00011) + |10001» + 
V5 

a 2 3 -= (|10100) + |01010) + |00101) + |10010) + |01001» , (9) 
V5 

where the normalization reads |ai,4| 2 + |a2,3 1 2 = 1- The task is then to find the optimal coefficients a M , such that C 
is maximized. 



B. Reflection symmetry 



One main big difficulty in the maximization problem of the concurrence is the growing number of free parameters 
a M for higher JV and p |l6j | . In this subsection we will elaborate on a symmetry argument that eliminates nearly half 
of the parameters. 

Recall the state under consideration, Eq. and the definition of the p-dimensional index /i, which specifies the 
number of steps from one up-spin to the next. It is clear that a cyclic permutation of the entries of /i will always 
result in the same state |/z), whereas an anti-cyclic one will result in a state where all spins are reflected around a 
certain site. Since is translationally invariant, it does not matter around which site the reflection is performed, 
and for convenience we can assume that the reflection interchanges the two sites that we are focusing on. Denoting 
the reflection operator by 7£, it is then obvious that g — > IZglZ leads to p' 2 ' — > [g^] (in our chosen basis) for any 
translationally invariant g. In terms of the matrix elements, reflection causes v,w,y — > v,w,y while z — * z*. In 
particular the reflected state has the same concurrence as the original one, C(lZglZ) — C(g). 
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Now, let us study the entanglement properties of g = (g + lZglZ)/2. We arrive at 



C(g) = \z + z*\- 2^-(v + v)-(y + y) = C(g), (10) 

and thus it is sufficient to consider states of the form g. Note that g commutes with the reflection operator, i.e. 
[g, 1Z] = 0, and therefore can be decomposed into projectors onto eigenstates of the reflection operator. By the same 
argument that allowed us to reduce the original problem to a pure state problem for each translational eigenvalue, we 
can now restrict our attention to pure states with translational eigenvalue 1 and odd or even reflection symmetry. In 
fact, using the observation from |jj that one only needs to consider positive coefficients a M in the optimal state, we 
immediately find that the even states will contain an optimal state. This result is non-trivial, intuitive as it may seem: 
A non-linear optimization problem could lead to solutions that break the underlying symmetries. Here, however, this 
is not the case and this fact reduces the dimensionality of the relevant parameter space significantly. 



III. ANALYTICAL RESULTS 



In the following we will consider various occupation numbers, starting with the trivial case p = 1: here the state is 
completely fixed and reads 

N-l 

l^(7v, P = i)> = -=^r i |^) , (ii) 

which is a generalized W-state. The concurrence between nearest neighbors (in fact, between any two sites) is given 
as C{N,p = 1) = 2/N. 



A. Maximal nearest neighbor concurrence for occupation number p — 2 

The case p = 2 is already much more involved. The concurrence of a given state can still be written in a compact 
way for all N, but it will turn out that the optimization problem can only be solved analytically for some small 
numbers of sites. The most general translationally invariant state for N sites with occupation number p — 2 is given 
by 

N-l 

\i>(N,p = 2))=Y. a ^Y. Tl \^) < ( 12 ) 

where, as explained above, fi is a index with two entries which describes the location of the up-spins in the ring; the 
first value for fi is 1, N — 1. As mentioned above, for the optimization problem we can take the coefficients to be 
real and positive. 

At this point we have to distinguish the cases of even and odd N: 

(1) Even N: 

The index fi stands for all inequivalent decompositions of N into a sum of two terms, starting at 1, N — 1, 
then 2,N — 2 and so forth until N/2, N/2. The last term of the superposition in Eq. I|12fl is special, because a 
translation by N/2 sites already leads to T^r \4>n/2,n/2) — \4 > N/2,N/'i)- Therefore this contribution has a different 
normalization, as mentioned below Eq. ||SJ. 

This becomes immediately clear from looking at an example; let us choose N — 4 and p = 2: 
\^{N = A,p = 2)) = oi, 3 |l,3)+02, 2 |2,2) 

N-l 

At=l,3;2,2 j=0 

= oi,3- Qiioo) + |0110) + |0011) + |1001>) + 

a 2 , 2 ^ (|1010) + |0101)) . (13) 



5 



For better clarity we will from now on simplify the above notation of the composite index jj, by dropping the 
second entry, i.e. di at-i = a>i and 1 1, TV — 1) = |1) and so forth. (This short-hand notation is of course only 
possible for the case p = 2.) Calculating the matrix elements for even N, one finds 

V = j^al , (14) 



2 

N 



n 2 
C = — max 

N 



(15) 



i (AT — 4 + af) for JV > 3 
otherwise 

N/2-1 

^2 a i a i+i + V2a N/2 -ia N /2 ] , (16) 

i=l 



The concurrence is thus given by the expression 

2\a 1 a 2 + a 2 a 3 H h a/v/2-20/v/2-i + V^a/v^-ia^l - \J N - 4 + o?, \ . (17) 



(2) OddN: 

Here the double index /i runs from 1, JV — 1 until (N — l)/2, (iV + l)/2. The normalization of all terms in the 
superposition is identical. An example of the most general state has already been shown for N = 5 in Eq. j9fr. 
With the short-hand notation introduced above we find 

y = ^al , (18) 

£(JV-4 + a?) for 7V> 3 
otherwise ' 

2 f^- 1 ^ 2 - 1 1 \ 

Z = N I a ' iai+1 + 2 a ( w -!)/ 2 I ( 20 ) 

The only differences to the case of even N are the last term in z being 2 " a (A'-i)/2 ms tead of V / 2«Ar/2-i a JV/2: an( i 
the upper summation bounds for z. The concurrence for odd N reads 

C = ^- max<2\aia 2 + a 2 a 3 H h a {N _ l) / 2 -ia{N-\)/2 + ^a(w-i)/ 2 l - \ai\^J N - A + aj,o\ . (21) 



From Eqs. (|17fl and (|21|l we see immediately that the first amplitude a% has a special role: it is the only one that 
appears explicitly in the last term of the concurrence, namely the term —y/vy. Note that the coefficient a\ specifies 
the only case where two neighboring spins are up. Considering both terms of the concurrence where a\ appears, we 
can find a simple argument when the coefficient a\ has to be zero in order to maximize the concurrence: for N > 8, 
the inequality 

2a 2 - yV-4 + a? < (22) 

holds for any a 2 and ai. Therefore, the concurrence is always increased by setting a\ — 0. Thus, for more than 8 sites 
(for both even and odd N) and the occupation number p = 2 one finds indeed the highest concurrence by using the 
constraint "no neighboring spins are up", i.e. C max is given by Eq. QJ. Note, however, that the actual coefficients of 
the optimal state are not trivial to find from the method in Q • Thus we have also used an alternative method which 
employs Lagrange multipliers. Our method does not only provide the value of the maximal concurrence, but also the 
coefficients a M that define the corresponding state. The details are described in the Appendix^ Using this method we 
find e.g. for the case N = 8 and p = 2 the results C max (N = 8,p = 2) = v3/4, where the coefficients of the optimally 
entangled state are given by a 2 = &3 = -\/l/2 and 04 = ^/l/3. 

For N < 8 we have in general a\ ^ 0, and the simple method from the Appendix cannot be applied. One has to 
maximize the concurrence explicitly. This was done for the cases A^ = 4 (here it turns out that one can reach the 
same concurrence as for p = 1, by taking a\ = \/l/3 and a 2 — -y/2/3); for N = 5 we have only two terms, namely 
\N = 5,p = 2) = ai|l,4) +<22|2,2) and arrive at the optimal solution a% — 0.298 and a 2 = 0.955. The case = 6 can 
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again be solved analytically: if the state is written as \N = 6,p = 2) = a\\l, 5) + a,2\2, 4) + a3|3, 3), the concurrence 
reads 

C(N = 6,p = 2) - | f^l-a\-al{ ai + V2a 3 ) - ^^2 + a^j , (23) 

where we have eliminated a 2 via the normalization constraint. The maximum can be found by demanding that the 
derivatives with respect to a\ and 0,3 have to be zero. This leads to the solution ax = and a,2 — C13 = \/l/2 for the 
maximal concurrence. Finally for N — 7 we have \N = 7,p = 2) = oi|l, 6) + d2|2, 5) + <23|3, 4) and 

C(N = 7,p= 2) = H ^2a 2 y/l- a\- a\- {a x - a 2 ) 2 + l- ai ^a\ + ^j , (24) 

where this time we have applied the normalization to eliminate 0,3. Calculating dC/dai, we find that a\ > ai => 
dC/dai < 0, whereas dC/da2 = is equivalent to a\ — 0,2 = (a| — a i)/ a 3- The last equality can only be fulfilled for (I) 
o>i > 0-2 > 03, where as noted dC/da < 0, for (II) a\ = a 2 = 0,3 which is clearly not optimal, or for (III) 0,3 > 0,2 > oj.. 

For the last case, we can use again an argument like for N > 8 to show that a\ — 0, 02 = v/(5 — \f§) /10, and 

a 3 = ^(5 + V5)/10 is optimal [13 and gives C max (7V = 7,p = 2) = (1 + \/5)/7. 

This completes the study of the case p = 2. To summarize, we have found a simple argument to show the solution 
in to be optimal for N > 8 and p = 2, and have explicitly performed the optimization for rings with seven or less 
spins. In the case N = 5, the seemingly simple one-parameter optimization results in a fourth-order equation that 
had to be solved numerically. All other cases were solved analytically. Note that N = 5 is the only non-trivial case 
where the optimal state contains two neighboring up-spins. We summarize the results for occupation number p = 2 
in Tab. □ (see also Fig. [TJ. 
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TABLE I: Maximal nearest-neighbor concurrence and coefficients for the optimal state, with p = 2. Only non-vanishing 
coefficients are given explicitly. 



B. Maximal nearest neighbor concurrence for 
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The case of half- filling, p = N/2, is special because the non- linear form of the concurrence © then is reduced to 

C{N,p = N/2) = 2max{|z| - |w|,0}, (25) 
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since here the matrix elements v and y are equal. As mentioned above, we may restrict ourselves to real and positive 
z, i.e. 2\z\ — z + z*. When this replacement and the equality \v\ — v is inserted in (|25[) . the concurrence can be 
written as the expectation value of a Hamiltonian. We start by expressing z = (01|g|10) and v — (00|p|00) in terms 
of Pauli matrices, 

z = -tr[g((7 + ®0] > (26) 
v = ~tr [g(a z - 1) ® (a z - 1)] . (27) 

Here cr ± = (o x ±i<r y ) denote the raising and lowering operators. Due to translational invariance, the matrix elements 
z and v are identical for any two neighboring sites. Thus, for non-zero concurrence and real z we arrive at 

C(N,p = N/2) = -tv(gH) , (28) 

where the corresponding Hamiltonian H is given by 

1 N 

H =w^ [~ <T?C7f + 1 " ° laV ^ + al ° 1 ^ + 1] • (29) 

i=l 

Here we have used that the magnetization per site is (of) = in this special case of half- filling. The Hamiltonian 1291) 
describes the Heisenberg model with nearest-neighbor interaction in an external magnetic field. (The relative sign 
between the z-coupling and the coupling in x- and y-direction does not change the physical properties of the system, 
as it can be absorbed by rotating the coordinate system of every second site by tt around the z-axis.) Note that the 
identity in H corresponds to an overall energy shift of 1/2. The ground state energy for the Heisenberg Hamiltonian 
in the limit N — ► oo is well-known 0] an d leads to 

C(7V_> oo,p = N/2) = 2 In 2 - 1 = 0.386 (30) 

As explained below, even for the case TV = 20 the concurrence is already very close to this value. 



IV. IS THE SOLUTION WITHOUT NEIGHBORING UP-SPINS A LOCAL MAXIMUM? 

The analytical solution for the maximal concurrence for any N and p was found in jjj by making the further 
restriction (ii) of no neighboring up-spins. This leads to y = in Eq. and one only needs to consider z. The 
matrix element z can be written as a quadratic form in the coefficients of Eq. JfjJ: z = & T Z&, where a denotes 
the vector with entries a M , and Z is a matrix, the dimension of which depends on N and p. Thus, the problem of 
maximizing the concurrence corresponds to finding the maximal eigenvalue of Z, i.e. 

C max (7V,p) = 2max{A fc (Z)}, (31) 

k 

where Ak(A) denotes the fcth eigenvalue of A. This diagonalization problem can be solved by the Jordan- Wigner 
transformation [l4| . 

It is in general a very hard problem to prove whether the solution in — we will call it OW solution from now 
on — is indeed optimal. In Sec. IIII Al we showed that for p = 2 and N > 6 no better solution can be found, but for 
general N and p even numerical approaches are very demanding due to the many free parameters and the non-linear 
character of theproblem. As we explain in this section, it is, however, much easier to numerically verify whether 
the solution in Q is at least a local maximum of the concurrence when the restriction of no neighboring up-spins is 
removed. 

What we basically need to do, is to expand Eq. J3) around the solution in |7J. In principle, the restriction (ii) leads 
to y — and thus to ^pvy — 0, which means that C as such is not a differentiable function of the coefficients of the 
state. Thus, one has to be slightly more subtle than simply finding the gradient of C. First of all, we observe that 
there is a natural grouping of the coefficients according to how many up- up pairs the corresponding state contains: 
the maximal number is p — 1, the minimal 0. A useful notation is thus to split the index on the coefficients a M in a 
"group" index n counting the number of up-up pairs and a second index j distinguishing the members of each group: 



(32) 
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The last notation is meant to emphasize that each can be viewed as a vector in its own right. The constraint (ii) 
allows only the coefficients in the 0-group to be non- vanishing, i.e. leads to a^ = (0, 0, ... , 0) T for tt > 0. 

If we restrict ourselves to cases where gcd(p, N) — 1, all normalization factors Mp in Eq. (JHJ) are equal to 1/y/N, 
and a given coefficient's contribution to y only depends on its group: 



0) = ]v 5>W 



(33) 



Since the number of up-up pairs gives directly the number of down-down pairs, a similar relation holds for v. 



As for z, things are naturally more complicated: z represents the mean value of 1 10) (01 1 and this flip of a down-up 
pair can at most create (or destroy) one up-up pair, i.e. there are cross-terms between tt and tt ± 1. In a natural 
matrix notation we get: 



z(a) = a T Za 



a n &i 



x P -i 



Zoo Zoi 

Zio Zn Z12 

Z21 Z22 





Zp~2 p-1 
Zp-l p—2 Zp-i p _x 



ao 
ai 



a p -i 



(35) 



(36) 



where the Z vn i are still matrices with indices running inside the groups. If we take the derivative of z at OW's solution 
a = w, we get 



dz 
dew. 



(w) = 2[Zw} 7 



2Stto [Zoow ]j + 25 n i [Z w -w ]j 
26-xo z(w) [w ]j + 25x1 [Z w wo]j, 



(37) 



where have used that Wo is an eigenvector (belonging to the largest eigenvalue) of Zoo- Since we have the further 
constraint |a| 2 = 1, only the component of the gradient perpendicular to the OW solution is interesting, i.e. only the 
second term in Eq. (|37|l needs to be considered. From the d„i we see that any local first order improvement of OW's 
solution can only be found by adding a small component from the space with exactly one up-up pair. Of course, when 
we do that y will no longer be zero, the square root in Eq. @ will grow and this will possibly cancel the benefits we 
could get from z alone. A priori we should test all directions in the 7r — 1 space, but since y does not discriminate 
between them, the one suggested by z suffices: if C does not grow when adding an infinitesimal component along 
•^io w Oj the solution of OW is a local maximum. The exact condition for this to happen reads: 



2|Zi W | < y/v(w) 



N 



(38) 



If, on the other hand, this inequality is fulfilled in the other direction, we know that OW's solution can be improved 
even locally and is thus for sure not the global maximum. 

In Fig. @ we will present the results of a numerical investigation of Eq. I|38|) for various N and p, together with 
the results of the non-perturbative investigation described in the following section. 



V. MAPPING TO A SERIES OF HAMILTONIAN PROBLEMS 



In Sec. lIII Bl we saw how the non-linear problem of optimizing the concurrence could be stated as a linear one in the 
special case of p = N/2. A method which exploits this idea further has recently been introduced by Wolf et al. [Tll |: 
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The maximal value of Eq. on translationally invariant states (for non-zero concurrence and z real and positive) 
coincides with the maximum of 

C(s)=z + z* - (sy+^v) , (39) 

where the maximization is done over all translationally invariant states and all s > 0. To see that this holds, simply 
note that the minimal value of sy + v/s with respect to s is 2^/yv, which is reached for s = vjy. Thus, after 
performing first the maximization of Eq. (|39[l with respect to s, one returns to the original problem. 

If on the other hand we optimize in the opposite order, then for a given s we can generalize the method from 
Sec. IIII Bl and relate C(s) to the expectation value of a Hamiltonian. Explicitly we have 

C(s)=-tv[ Q H( a )}, (40) 

where 

1 N 

i=i 

This Hamiltonian corresponds to the quantum XXZ-model with an external magnetic field B(s) := s — 1/s and an 
anisotropy parameter A(s) := (s + l/s)/2. For s = 1 Eq. H41[) reduces to Eq. H29|) . The XXZ-model is completely 
integra ble. and in principle all eigenvalues and eigenstates of the Hamiltonian l|41l) can be found using the Bcthc 
ansatz [l5|. Unfortunately for general p, corresponding to general magnetization, the solution is very difficult to 
handle even in the limit N — > oo [l^, an d we have to resort to numerical solutions. 

For a given s, we minimize Eq. I|41|) by direct diagonalization. It turns out that the dependence of the ground-state 
energy on s is rather simple, such that the final optimization with respect to s provides only a modest complication. 
In Fig.^we plot the results obtained for all relevant p as N ranges from 1 to 24. "Relevant" p are those with p < N/2, 
as the concurrence is identical for occupation number p and N — p. For a fixed p the maximal concurrence has in 
general the following behavior as a function of N: from the first relevant N — 2p, it slightly decreases and then grows 
to a maximum around N = 3p, after which it tends slowly to as N — > oo. 




FIG. 1: Maximal nearest-neighbor concurrence for N = 4, . . . , 24. Each line and symbol type corresponds to a fixed p: the 
left-most curve illustrates the case p = 1 and is marked by the symbol "+". The second curve from the left shows p = 2, and 
so forth. The explicit legend is also given in Fig. [3] Curves with higher p's peak at higher TV, the optimal value scaling roughly 
as 3p. 

In order to analyze more carefully when OW's result is optimal, in Fig. [2] we show parts of the C max (iV) curves 
for the OW constraint of no neighboring up-spins, for the improvement via perturbation theory, and for the exact 
solution. As can be clearly seen, optimal states with no neighboring up-spins can be found for small numbers of sites, 
namely for N smaller than the value for which the concurrence has its maximum for a given p. Note also that the 
local optimization via perturbation theory and the exact calculation always start to deviate from the OW curve at 
the same N. This implies that when the OW state is locally optimal, it is also globally optimal. 
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N 



FIG. 2: Maximal nearest-neighbor concurrence for p = 2, . . . , 7 as function of the number of sites: OW result (crosses), 
perturbation theory (circles), and exact solution (squares). Note that the perturbation theory and the exact solution in all 
cases meet the OW result at the same N. This implies that OW's solution is never a local maximum without being also the 
global maximum. 

It is furthermore interesting to study the values for the optimal s as a function of N and p, see Fig- EI These values 
immediately give the physical interaction that leads to the ground state with maximal nearest-neighbor concurrence. 
Notice that s — 1 for p = N/2, and that for any value of p the optimal value of s goes to infinity with N — > oo. This 
indicates that the OW solution is indeed optimal for large N, and that the occurrence of neighboring up-spins for the 
maximization of the concurrence is a true finite size effect. 




1 1 1 1 ' — 

5 10 15 20 

N 



FIG. 3: Optimal value for the parameter s as a function of N and p. 

Let us also mention the exceptional case N — 8. In 0, no solution with a concurrence above the asymptotic value 
C max (N — > oo) w 0.434 was found for any p. Is this an artifact of the restriction to no neighboring up-spins, or a 
true finite size effect? Our studies show that the latter is true: we find C max (iV = 8,p = 3) w 0.431. Thus, together 
with the results presented in Table [I] the occupation number p — 2 is shown to be optimal for N = 8. However, 
C max (7V = 8,p = 2) = \/3/4 w 0.433 < C max (iV -> oo) holds — a peculiar feature of a spin chain with 8 sites. 

VI. BEYOND NEAREST NEIGHBOR ENTANGLEMENT 

A natural extension to the analysis done in the previous section is to examine entanglement between spins separated 
by at least one site. In principle, the formalism of Sec. IIII Al can be applied to the maximization problem of the 
concurrence between spins separated by any number of sites. This approach yields results similar to those of Sec. IIII Al 
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For example, when maximizing the next-nearest neighbor entanglement for p = 2, we find that using the constraint 
"no next-nearest neighbors are both up", i.e. taking a 2 = 0, does not decrease the concurrence for N > 8. 

In this section we will treat the issue of reducing the gth-neighbor maximization problem of a ring with N sites 
to a nearest-neighbor maximization problem of a ring with less than N sites. In contrast to the previous sections, 
we will not fix the number of up-spins; we simply look for the maximum of C with translational invariance and no 
superpositions of different p. 

Assume first that q divides N, i.e. q\N. From now on we will denote the maximal gth-neighbor concurrence of a 
ring with N sites by C lmax ( N, q) . We denote by liV) 013 * a state that reaches this value and by p opt the corresponding 
number of up spins. 




FIG. 4: Decomposition of the ring into smaller rings, for the case N = 8 and next-nearest neighbors. 



Consider now the decomposition of the large ring into q smaller rings. Each of them consists of every qth qubit 
of the large ring and has N/q sites. The example of a ring with 8 sites and next-nearest-neighbor entanglement, i.e. 
the case N = 8 and q = 2, is illustrated in Fig. 0] The state that maximizes the nearest-neighbor concurrence of a 
small ring, C max (N/q, 1), will be denoted by \N/q) opt . The q-fold tensor product of this state is a candidate for the 
state |7V) opt : it is translationally invariant since T can be decomposed into a permutation of the set of small rings 
and translations within them. We thus have: 

C(|A7g) opt « • • • ® \N/q)°v\ q) < C maK (N, q) = C(\N) opt ,q) (42) 

" v ' 

q times 

To show than one cannot do better than using this candidate, we note that the entanglement between the first and 
the qth. site of the state |7V) opt is calculated by tracing over all states but the first and gth. In particular, we trace 
over all but one small ring. Let us denote the state obtained by tracing out the q — 1 irrelevant rings by q^ SmT '\ i.e. 

=tr othcrrmgs |A0(AT pt ■ (43) 

Therefore, 

C{\Nr\q)=C{e^\l). (44) 

Now, because the full g was translationally invariant, g( s ' r ') is also translationally invariant. This means by definition 
that 

C(g (s ' r - } , 1) < C max {N/q, 1) = C(\N/q) op \ 1). (45) 

Since the tensor product of Eq. I|42|l saturates this inequality, it must be optimal. We conclude that when q \N, 

C max (A^, q) = C m&x (N/q, 1). (46) 

Let us now drop our assumption q\N and assume that N and q are coprime, gcd(iV, q) = 1. If one then tries to 
construct "small rings" by taking steps of q sites, one will eventually visit all sites, i.e. one only finds one such ring. 
This simply reflects the fact that T q generates the same group of translations as T does. Then it is clear that by 
reordering of the spins one can construct a ring on which all spins that were formerly separated by q sites are now 
adjacent and on which translational invariance is equivalent to translational invariance with respect to the original 
ordering. Therefore, when gcd(7V, q) = 1, 



C max (7V,g) = C max (AT,l). 



(47) 
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In general, constructing small rings by taking steps of q sites will lead to gcd(7V, q) rings, each of length N = 
N I gcd(iV, q). On these rings the former qth neighbors are now nearest neighbors, and an optimal state \N, q) opt can 
be constructed similarly to the q \N case above as a tensor product of small ring states \N, l) opt . In this way, the 
problem can always be reduced to the case of nearest neighbors with the result that for any N and q: 

C max (N,q) = C max (iV/gcd(7V, 9 ),l). (48) 

The method of this section naturally generalizes to other symmetry groups acting on N qubits: If we choose to 
optimize a given pair's entanglement, the symmetry group will give some restrictions, but these may "factor" so that 
tensor product solutions of smaller problems can be applied. In the ring symmetry case, Eq. I|48|) shows that the 
reduction is maximal in the sense that there is essentially only one class of problems, namely nearest neighbors on 
a ring with translational symmetry. This is due to the particularly simple structure of the group of translations; in 
general one symmetry group may give rise to many fundamentally different optimization problems. 



VII. DISCUSSION AND SUMMARY 



In this paper we have studied maximal bipartite entanglement in translationally invariant spin chains with periodic 
boundary conditions, where our only restriction was to consider eigenstates of the z component of the total spin. We 
have concentrated on bipartite entanglement, as entanglement in this case is fully understood, contrary to multipartite 
entanglement. 

Naturally, some related questions arise: 

• What is the maximal entanglement between one spin and the rest of the ring? 

• What is the possible structure of three- or multi-party entanglement in the ring? 

• What is the optimal nearest-neighbor entanglement on a 2-dimensional lattice? 

Here we will only answer the first question: states with fixed occupation number p are the eigenstates of the z 
component of the total spin, i.e. they are invariant under a global rotation around the z-axis, U = ®i!=i °f • This 
leads to [a z , g^} = 0, where g^ denotes the one-particle reduced density matrix, which is identical for all sites due 
to translational invariance. Thus, g^ is diagonal, i.e. g^ = (1 + s z a z )/2. Here s z = (er z ) is the magnetization in 
^-direction for one site. The explicit form of the reduced density matrix is therefore 



y N 



(N-p) 
p 



(49) 



Note that this form is independent of the choice of coefficients a M , and results only from translational invariance and 
the condition of fixed p. The von Neumann entropy of this density matrix is a good measure for the entanglement of 
one qubit with the rest of the chain. Therefore, we find maximal entanglement in this sense when p = N/2. 

In summary, we have investigated the maximal nearest-neighbor concurrence in translationally invariant entangled 
rings, where the number of up-spins p is fixed. We presented analytical results for small N and p = 2. It also turned 
out that the solutions in can be improved by simple perturbation theory. We then used a linearized version of the 
concurrence, which simplified the numerical calculations, and with this method determined the maximal entanglement 
for any p with N < 24. The maximal entanglement was shown to be related to the ground state energy of an XXZ 
Hamiltonian. Let us emphasize our observation that a state with maximal nearest-neighbor entanglement corresponds 
to a Hamiltonian with only nearest-neighbor interactions. For small numbers of sites N and fixed p we gave the explicit 
structure of the Hamiltonian which leads to the state with maximal entanglement of nearest neighbors. We found 
the peculiar result that for the exceptional case N = 8 the maximal entanglement is below the entanglement in the 
thermodynamic limit. Finally, we reduced the problem of finding the maximal entanglement between spins which are 
separated by more than one site to the problem of nearest-neighbor entanglement. 

This work has been supported by the Deutsche Forschungsgemeinschaft (SFB 407 and Schwerpunkt "Quanteninfor- 
mationsverarbeitung" ) , the European Union RTN "Cold Quantum Gases" , and the Danish Natural Science Research 
Council. 



APPENDIX 



Here we describe a method to find the maximal concurrence and the according state for the case N > 8 and p = 2. 
We take the normalization constraint, namely af = 1, into account via Lagrange multipliers. 
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For even N we then have to solve the system of equations 

a 3 + 2a 2 A = , 
a 2 + a 4 + 2a 3 A = , 

&3 + &5 + 2<24A = , 

&N/2-2 + V2apf/2 + 2a N /2-\\ = , 

V2a N/2 -i + 2a N/2 X = 0. (A.l) 

where A is the Lagrange multiplier. For odd N the last two equations in the system of equations look slightly different: 

+ a 3 + 2a 2 A = , 
C12 + 0,4 + 2a 3 A = , 
a 3 + a 5 + 2a 4 A = , 

a (JV-5)/2 + A(iV-l)/2 + 2a( A r_ 3 )/ 2 A = 0, 

a (JV-3)/2 + a(iv-i)/2 + 2a( A r_ 1 )/ 2 A = 0. (A. 2) 

In both cases we now multiply the first equation with et 2 , the second with 03 and so on, and then add all equations. 
This leads to the expression 

C,„„ . -« (A.3) 

for the maximal concurrence. Note that this equality holds for both even and odd N. By inserting the first equation 
in the system of equations i|A.l|) into the second, the new one into the third, and consecutively until the last equation, 
one finds a polynomial equation for A of order N/2 — 1 for even TV. By proceeding with (|A.2|I in the same way one 
arrives at a polynomial equation of order (N — 3)/2 for odd N. For small N, these polynomial equations can be solved 
analytically, and for larger N one can still find the zeros of the polynomial numerically. 

Let us illustrate this method with a simple explicit example, namely the case N = 8 and p = 2. This is also 
a particularly interesting case, because in Q in was shown that the maximal concurrence for this case (with the 
constraint of no two neighboring spins being up) is smaller than for the thermodynamical limit N — > 00. 

The system of equations for N = 8 and p — 2 is 

a 3 + 2a 2 A = , 
a 2 + v2a4 + 2a 3 A = , 

V2a 3 + 2a 4 A = . (A.4) 

By inserting these equations successively into each other we arrive at 

4A 3 - 3A = . (A.5) 

In order to maximize the concurrence in Eq. (|A.3|I we have to find the minimal solution for A, which is A = — \/3/2. 
Therefore C(N — 8,p = 2) = \/3/4. The coefficients of the optimal entangled state are then easily found to be 
a 2 = \A/6; a 3 = \A/2 arL d (14 = \/l/3. The results of other cases of small N are given in the main text. 
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